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The effective QCD Hamiltonian is constructed with the help of the background perturbation theory, 
and relativistic Feynman-Schwinger path integrals for Green's functions. The resulting spectrum 
displays mass gaps of the order of one GeV, when additional valence gluon is added to the bound 
state. 



^ ' Mixing between meson, hybrid, and glueball states is defined in two ways: through generalized 

, Green's functions and via modified Feynman diagram technic giving similar answers. Results for 



mixing matrix elements are numerically not large (around 0.1 GeV) and agree with earlier analytic 
estimates and lattice simulations. 



o 

Cu'. 1 Introduction 

^ ; Gluonic fields play a double role in the dynamics of QCD. On one hand they 
create the QCD string between color charges which substantiates confinement. 
^ ■ On another hand valence gluons act as color sources and can be considered 
c3 ; as constituent pieces of hadrons on the same grounds as quarks. This is a 
QCD string picture of hybrids and glueballs developed analytically in and 
confirmed by recent lattice calculations [0] . On experimental side the study of 
hybrid and glueball states is not yet conclusive, and more experimental 
and theoretical work is needed, especially a quantitative treatment of mixing 
between mesons, hybrids, and glueballs. 

One typical feature of the QCD string approach, supported by lattice data 
, is that the addition of every valence gluon in a hadron increases the hadron 
mass by 0.8 - 1 GeV, and this is true also for purely gluonic states - glue- 
balls. Thus a large mass gap exists between meson ground state and its gluonic 
excitation, which makes it possible to consider gluonic admixture in a given 
hadron state as a perturbation. Therefore in the zeroth approximation one has 
a Hamiltonian for the diagonal states of the fixed number of quarks q and 
valence gluons while in the next approximation one calculates the mixing 
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between the states perturbatively (unless masses of the states happen to be 
almost degenerate, in which case one solves a matrix Hamiltonian). 

The mixing between meson and hybrid states was considered previously an- 
alytically in In []5| the hybrid wave function was taken in the cluster 
approximation, with gluon wave function of bag-model type factorized using 
the qq potential model wave function. The resulting mixing matrix elements 
(MME) are not large, supporting the iterative scheme described above. A sim- 
ilar approach, based on the nonrelativistic constituent quark and gluon model, 
was used in [g] to study mixing between l^ccg hybrid meson and charmonium. 

The mixing between glueballs and other hadrons was studied on the lattice 
[0], with reasonably moderate MME. 

It is a purpose of the present paper to develop a general formalism of the 
QCD bound states including valence gluons, with the special attention to the 
mixing between states. This formalism is based on the QCD string approach 

and applicable to both relativistic and nonrelativistic systems, massive or 
massless quarks and gluons|^. 

Two features are important for this formalism. First, it is derived directly 
from QCD with few assumptions, supported and checked by lattice computa- 
tions. Second, the only input of the theory is the fundamental string tension cr/, 
current quark masses (renormalized at the scale of 1 GeV), and strong couphng 
constant ag. (For the total hadron masses one needs to subtract a selfenergy 
for each quark and antiquark equal to 0.25 GeV, but this does not affect wave 
functions and mixings which will be the main goal of this paper). 

The plan of the paper is as follows. In the next section we develop the 
background perturbation theory to separate valence gluons from the background 
gluonic fields forming the QCD strings. On this basis the diagonal part of 
the QCD Hamiltonian is written for quarks, antiquarks, and valence gluons 
interacting via string connection and main characteristics of the spectrum are 
established for mesons, hybrids, and glueballs. 

In section 3 the part of the Hamiltonian responsible for the mixing of mesons, 
hybrids, and glueballs is identified and MME is written in terms of solutions of 
the diagonal Hamiltonian, using formalism of modified plane waves. 

In section 4 concrete calculations of MME are presented for meson-hybrid 
case based on the Green's function formalism. An analogous treatment of the 
meson-glueball case is given in section 5. 

short version of an approach in the same direction with a different form of matrix elements appeared 
recently ^. Results of (|] are similar to ours qualitatively, with some numerical differences. 
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In concluding section discussion is given of the results obtained in this paper 
and comparison is made with calculations with ohter model, lattice simulations, 
and experimental data. 



2 Green's functions and Hamiltonian 

We follow in this section the procedure developed in detail in [0, |TD[ ( for earlier 
refs. see [|rT|) and therefore here we recapitulate only the main points. The total 
gluonic field is represented as a sum 

A^ = B^^ a^, (1) 

where 5^ is the nonperturbative part and - perturbative (or valence gluon) 
part, which shall be treated in form of the perturbation series in ga^. One can 
formally avoid the problem of doublecounting using the 't Hooft identity, 
which allows to represent the partition function as 

Z = ^J e'^^^^DqDqDA = ^J r]{B)e~^^^+''^ DBDaDqDq, (2) 

where A^^, A^^' are normalization constants. 

Here S is total Euclidean action and t]{B) - an arbitrary weight of averaging 
over DB. Expanding S{B + a) one obtains 

S{B + a) = S{B) + Si{a, B) + 53(0, B) + 53(0, B) + S^{a) + Si{a, q, q) (3) 

In (0) Sn{a, B) denotes terms of power n in a^, and Si{a,q,q) is the mixing 
term, which will be of main interest in section 3. The terms and <5'4 contain 
additional powers of g and can be treated perturbatively, the term Si{a, B) was 
considered in [|rD| and shown to yield a small correction to the leading terms 
S{B) and 5*2(0,5), the latter defining the valence gluon propagating in the 
background B^ with the Green's function (in the background Feynman gauge 

H-H) 

Gf. = (^aV + '2i9F,.)-\ D, = dx- igB,, f;^ = d,Al - d,A; + gr'^'A'^Al. 

(4) 

The explicit form of all terms in (^ and of the corresponding Green's functions 
is given in [j^, |T3], here we only quote the results for the Green's function of 



the state containing quark, antiquark, and any number of valence gluons in the 
leading A^^ approximation: 

n 
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For (multi)glueball Green's function one has a similar representation with miss- 
ing factors Sq^Sq. Here Tin^Tf are initial and final vertex operators and {...)b 
means averaging over fields with the weight t]{B), and quark Green's func- 
tion Sq can be written as 

S, = -i{D + m)-' = -i{m-D){m^-Dy = -i{m-D){m''-bl-a^,F^,)-\ 

(6) 

In what follows we shall treat spin-dependent terms for gluon and quark Green's 
functions in (g) and (g) as a small correction, which is supported by exact calcu- 
lations and comparison with experiment and lattice data HI. Correspondingly 
one can use for the quadratic in parts the Feynman-Schwinger (world-line) 
path-integral representation (see [[12], [T^ and refs. therein) 



^ 2 
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D%1 = /f ds j (Dz),,e-'' exp(»g B,dz,), K = \ £(^?dT. (7) 



It is important that background field enters ([7j) only in the exponentials 
and one can show [jT2j that in the total Green's function (^ in the leading Nc 
approximation all these exponentials combine into products of Wilson loops. 
Here we make another assumption, supported by lattice data, that Wilson loop 
has a minimal area law behaviour with the string tension af (for mesons and 
hybrids) and (Jadj = |cr/ for 2g glueballs. Thus one reduces the system of g, g 
and n gluons to the open string with the ends at q and g, and n gluons "sitting" 
on the string, and glueballs reduce to an open adjoint string (or equivalently 
for large - a closed fundamental string). 

At this point one can define the Hamiltonian of the system H, using the 
relation for the general Green's function of the type 



GiX,Y) = {X\expi-HT)\Y) (8) 

where T = X4 — Y4. 

The Hamiltonian H can be defined on any hyperplane, in case of center-of- 
mass system one gets the Hamiltonian which was obtained for the qq system in 
and generahzed in [jT?! to hybrid and in []T^ to the glueball case (see also 



r5|j for earlier papers on hybrid Hamiltonian). 

Here we quote the simplest version of the Hamiltonian for qq and n gluons, 
where string rotation contribution to the moment of inertia is treated as a 
perturbation: 

H = Ho + AHl + AH, + AHa (9) 
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where we have defined. 

H, = !^+t^+ ^ Ei±Z!!k+^^{|r,_r,| + |ri-r,| + ... + |r„-r,-|}. 

^ fc=l ^ k=q,q,l...n "^l^k 

(10) 

Here irik = for k = 1, ...,n and fik defined as (see [14]) are to be found 
from the minimum of the Hq, which is approximated within 5% accuracy by 



the minimum of the eigenvalue M^'''^ [0, Hq"^^ = M^"^"^, 



AHl gives zero contribution when all interparticle angular momenta are zero 



=0' k = q,q, 1, n; z/ = 0, 1, 2, ... (11) 



and in particular for the meson case has the form [jT4 



The Coulombic part AHc takes into account lowest-order Coulomb exchanges 
between quark and antiquark. It is argued in [TBI that Coulomb exchanges 
between valence gluons are strongly suppressed by higher-loop corrections and 
hence can be neglected in the first approximation. 

Finally, the spin-dependent term AHg has the following form for qq or gg 
case 

• Sfc S, • L, 62 ^ g J . 3(Si • n)(S2 • n) - Si • S2 
AHs = a 2^ h6i2. 2 — ^ 2^S^-Lj+c . 

(13) 

In the c.m.s. Li = — L2 = L 0,61,62, and c depend on distance r between 
qq{gg) and can be found in 

Let us now discuss the properties of the resulting Hamiltonian (|9]). First 
of all, one should stress that it is a fully relativistic Hamiltonian. Indeed, 
neglecting for the moment the interaction term in (PUI ) and finding fik from ([TT|) 



one immediately obtains /ik = v"^fc + P ' ^^^^ ^^^^ f^k plays the role 

of the relativistic energy of a quark or gluon. Moreover, the form of AHs is 
not a result of expansion in inverse powers of large masses fii, fik, but rather is 
a result of Gaussian approximation for the average of spin-dependent factors 
{llnexp{gal]^^F-!!^))B (see [jlj] for details). 

It is important to stress that m^, rriq entering the Hamiltonian ([T^) denote 
quark current masses, renormalized at the scale around IGeV, and nowhere we 
use as input constituent masses of quarks or gluons. 
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The spectrum of the Hamiltonian (|^) was calculated for many systems, in- 
cluding mesons, hybrids, and glueballs, for a review see As a recent example 
see calculation of gluelumps in JT^ . 

The characteristic feature of the spectrum is that each constituent (quark or 
gluon) contributes to the total mass a quantity approximately equal to 
where fi^^^ depends on the number of eigenvalue u and is expressed through 
aj (or dadj for glueballs). nf'' changes from 0.35 GeV for massless quarks in 
the lowest meson states, and = 0.52 GeV for gluons in the glueball. 

Therefore any gluon in a hybrid adds around 1 GeV to the total mass, and the 
same is true for glueballs: the lowest three-gluon glueball is approximately 1 
GeV heavier than the lowest two-gluon glueball. This fact supports the idea 
outlined in the Introduction that the diagonal eigenvalues of the total Hamil- 
tonian have energy gaps around 1 GeV and it may be a good approximation 
to treat mixing due to valence gluon excitation as a perturbation provided the 
MME is much less than 1 GeV. 



3 Mixing matrix elements 

We turn in this section to the part of interaction, Si{a, q, q), which is responsible 
for the mixing between hadronic states, differing by the number of valence 
gluons. It has the form 

Si{a,q,q) = 9 J q{x)d{x)q{x)d'^x. (14) 

Here d{x) = a^t'^pjfj,: and qa, qp all have color indices, whereas in the Hamil- 
tonian ( pXl| ) and its eigenfunctions the color indices are absent because of color 
averaging in Cqq^^g^ (|^) and hence in H. 

To understand how the matrix element of Si is taken between colorless 
hadronic eigenfunctions of H ^ one can use the formalism of Green's functions, 
described in Sections 4,5, and corresponding to the diagrams in Figs. 1-4. 

Here we shall choose a simpler way, which leads to the same results as in 
sections 4,5, but more familiar for the reader accustomed to Feynman diagrams. 

The rules are simple: represent each q{x)^q{x) and a^(a:) by equivalent col- 
orless fields with familiar plane- wave expansion 

^(x, t) =J2 -f^=7T^i(^(^^ exp(zk • X - iiiqt)+u''{k, (j)dX exp(-zk • x + iiiqt))] 
k V2/^9(k)^ 

(15) 
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= E /o \ NT^ [^^P(^^ ■ ^ ~ ^/^^)eL'^^CA(k)+e{;^)4(k) exp(-zk • x + z/it)], 



(16) 

where 6a and c\ are annihilation operators for the quark and gluon respec- 
tively, dx is creation operator for antiquark and u{k) is normalized according 
to the condition u{k)u{k) = 2fiq{k), V - is the 3d volume and e^^^ is the gluon 
polarization vector. 

One can easily check at this point that introduction of ([T^) into the free part 

of the QCD Hamiltonian — ^ — immediately reproduces the gluonic part of the 
effective Hamiltonian ( PI3|) (the same is true for q if one uses the quadratic in 
k form of the Hamiltonian). 

In (|T5|), (|T^) /ig(k) and /i(k) are the same quantities as in (pXI|) , to be defined 
by the minimization procedure in the system, where quark or gluon enter as 
constituents. 

Another important point is the number of independent polarizations of a 
bound gluon. We shall use as in y-[lTT| "the gauge-invariant gauge condition" 

D^a^ = 0. (17) 

For a valence gluon in the gauge-invariant hybrid wave function one has 

*^(a:, z, y) = g(x)$(x, z)a^{z)^{z, y)q{y). (18) 

Here $ are parallel transporters depending on 5^, $(a:, y) = Pexpig Jy B^dz^. 

Differentiating ^ ^ in one obtains additional term due to differentiation 
of the end points in $(a:, z), $(2;, y), so that one has 

d 

^*^(x, z, y) = q{x)^{x, z){D^a^{z))<^{z, y)q{y) + ... (19) 

where dots imply the terms from the contour differentiation, i.e., due to addi- 
tional gluonic excitation of the string. 

In the diagonal approximation (and keeping in mind that gluonic excita- 
tion implies a gap of 1 GeV), we disregard those terms, and hence have the 
(approximate) condition 

^^(f)(x,z,^) = 0. (20) 

This means that only 3 gluon polarizations are physical, i.e., the same situation 
as for an off-shell photon or gluon. In what follows we shall retain only 
/i = 1, 2, 3. It is clear that the same reasoning applies to gluon in glueball. 
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We are now in the position to write the general form of wave functions for 
mesons, hybrids, and glueballs. In the momentum space and in the second- 
quantized form they can be written for a meson 

0^"'^ = fL::a7\pi. P2K{pi)dM. (21) 

where o'i,cr2 are polarizations (helicities) of q and and ai,a2 are Dirac 4- 
spinor indices. 

For a hybrid one has 

^^""^ = ft,a:f \Pl, P2, P,)btMKAP2)ct{P3) (22) 

where A is the gluon polarization and /i is discussed above, it enters as in e^^^ 
in (|TB|). Finally for the glueball one has 

= 4^,^^•'(Pl, P2)4,(pi)4.(P2). (23) 

One should note, that all wave functions ([2Tl) -( p3|) are given in the represen- 
tation when the total angular momentum J and its projection Mj are not 
projected out. In fact the operator (PIp is a 16 component structure, and one 
can use the classification introduced in to distinguish positive and negative 
energy states of quarks using the so-called /O-spin, p = ± and usual spin states 
for each quark. From spin and /o-spin states one can construct the states with 
given total momentun and parity []22|. These scheme was exploited and devel- 
oped in the series of papers where the qq interaction was used containing 
both scalar and vector confining parts. 

A similar scheme can be used for hybrids, but the resulting calculations 
are rather cumbersome. Therefore in this section we only list some schematic 
expressions with coefficients for mesons which can be found in []2^ and for 
hybrids not yet available (to the knowledge of the author). Thus one can 
write instead of wave functions (PT|)-([2^) the functions and ^j^^ 

with given total angular momentum. Each of this functions consists of several 
components differing in total spin, orbital momentum, and p-spins. It is actually 
those combinations which should be inserted everywhere below in this chapter 
instead of ^^^\ and ^'^^\ respectively, but we keep for simplicity reason 
the functions (pi|)-([231), referring the reader to the Appendix for more details 
and discussion. 

We are now in position to write down the MME between states ([2T|) -(|23|), 
using the interaction Hamiltonian, obtained from ([T^l), namely 



Hi= g I g(x, 0)a(x, 0)g(x, 0)rf^x = 
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= g E g-(k0<i{k3k{k,)^-|%^(2.)¥^)(ki + k, + k3). (24) 
k„k„k, (2'^) (^tt)-* (27r)-' 

For the meson-hybrid mixing one has 



The normahzation condition for $(^^),$(^) looks hke 

(27) 

where we have defined 

dP3 = n ^,{27Tf6^'\p, + P2 + Pb), (28) 

^^2=n^(27r)¥3)(Pi + P2). (29) 



.^1 (27r)3 

Going from momentum to coordinate representation and suppressing for the 
moment the momentum dependence of the factor r[pi,p[, k^) in ([251) , where 



r(pi,y„fc3) = ^=4=, (30) 



one has 

($(^^), //i$(^)) = gj dW\r)rip^''\0, r) (31) 
where we have defined 

^^''\r) = / /l:rHPi,P2)e^P^-'^+^P^''^^i^2 = 

= / fL?,a?\Pl. -^^y^'''-§^2^ r = - (32) 

^^""HO, r)=l ft,a?f Kpi,P2, P3)e^P^-"^+^P^-"^+^P3-r3^P3 (33) 

with ri — ra = 0, ri — r2 = r. 

In the next two sections we shall use another formalism to derive MME - 
the formalism of Green's functions, which enables us to use our Hamiltonian 
technic described in chapter 2. 



Figure 1: Meson Green's function describing propagation of the quark and antiquark (solid 
lines) from the point x to the point y. The hatched interior of the figure implies presence of 
nonperturbative fields in the form of the fundamental string world sheet. 

4 Mixing between meson and hybrid states 

Consider a nonsinglet qq state; the corresponding Green's function in the 
quenched (large Nc) approximation can be written, according to general for- 
mulas []T^, |r3|, as (see Fig. 1 for the corresponding Feynman diagram) 

G,,{x,y) = (v5(x)r(-*)v^(x)v5(2/)r(-)v^(2/))5,„^^^^ = 

= {trr^'^'''^S{x,y)r^^^S{y,x))B,a. (34) 

Here averaging over 5, a is assumed with the action given in (0)5(0) and S{x, y) 
is the quark Green's function, 

Six, y) = -i{D + m);^y, D^ = d^- igB^. (35) 

Writing S = —i{m — D){w? — D"^)'^, one can reduce the Wilson-loop path 
integral for Gqq to the form 

Gqq{x,y) = (tre«r(°-*)(m-I)) f f dsdse-""-^ {Dz),y{Dz),yr^'-\m-b)WF) . 

(36) 

Here Wp is the Wilson loop with insertion of operators aF, defined in and 
trace is over color (c) and Lorentz [a) indices. 

Our primary task now is to transform the integral / ds{Dz)xy of each quark 
Green's function as follows: 

d.(D.)., = d.ng|l^e'>(--E.^<"», Ne = s. (37) 
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Introducing now the parameter /i(z4) playing the role of constituent quark 
mass according to 

where 2fls = En Az4(n) = T,n ^s2fi{n), one obtains 

DsiDz)^, = ^^e'P-<^-y-S.. ^(""D^.Dm. (38) 

Thus the path integral for each Green's function acquires the 3-d form, which 
we call G(3(i), with all Zitterbewegung (z-graphs) contained in the integral over 
D/i. If one does, as we usually do for all systems (except pion and kaon), the 
stationary point procedure in integration over Dfi, then one finds the smooth 
"constituent mass trajectory" fio{t){t = 2:4), or a simple approximation to it, 
the constant constituent mass /io, found from the minimum of the Hamiltonian 
eigenvalue In this approximation (yielding accuracy around 5% for masses 



TB|] ) one can identify fi in ([38|) and fiQ. 



Thus one can write 

S{x,y) = ^^^G3,(x,2/) = ^^^(x|n)e-^^"l--^l(n|^), (39) 

where M„ is the n-th eigenvalue of the Hamiltonian. Similarly for the gluon 
Green's function one writes 

G,,{x,y) = - 2iF,Xl = ^{x, ^\n)e-''-\--y\{n\y,v). (40) 

The form in (P^) , (|iD|) is highly symbolic, since quarks and gluons do not 
propagate separately (and do not have separate eigenfunctions (|n)), but rather 
form the common bound state. In the free case one can identify Jlq with the 
energy Eq, Jig with tOg, and \n) = so that (|39|), (|40| ) go over into well- 

known representation 

c/ N ■/ f exp(ip-x-^|x4-^4|) ..^s 
5(x, ,) = -^{m-^)j ^ , (41) 

r d'kexp{i{k-^-u\x^-y,\)) 
G,.(x, y) - V J (27r)32a;(A:) ' ^^^^ 
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Consider now the hybrid Green's function 

G('')(,t("„t<^), yP), yW) = (*(")+(,T<'), x<^), x(=")*('')(j/<'), y<^), j/<^)))5,„,,. 

(43) 

Omitting parallel transporters for simplicity, one can write 

where we have denoted A^^'^^ = ^^^^ for quark and antiquark respectively. 

The Hamiltonian eigenfunctions (nly'^^^ y'^^^ y^^^) are characterized by the 
cm. momentum P, and boundstate quantum numbers n and the sum over n 
in ( ^) can be written as 



where $y^]r}^'^^ are Jacobi coordinates, while X, Y are cm. coordinates. 

We are now in the position to write the amplitude, corresponding to the 
Feynman diagram of Fig. 2, where the valence gluon is emitted at the point 
u and absorbed at the point v. The corresponding term in the Lagrangian is 
given in ([Tip , 

Si{a,q,q) = j q{x)ga{x)q{x)d'^x. (46) 
12 



The general form of the amphtude of Fig. 2 is 

X S{y, v)gj^t^d%S{v, x)) b- (47) 
Here trace is over both color (c) and Lorentz [a) indices. 



Introducing now the representations as in (^), one obtains finally 



2/i,(n)(Mr - Mm){Mm - Mi")' 



X 



X i?!5_e*P-(x-y)-i?M(y4-x4) 
(27r)3 

Here notations are used 



(48) 



V^^l = gj^M{vri^^{Q.v)dh (49) 

and v^M(r) , '^V^n are meson and hybrid eigenfunctions respectively (with the 
gluon in the latter having polarization /i), while Mm^m]^^ are masses of the 
meson and hybrid (in the n-th excited state) respectively. Note that while Gqq 
(|3§) is 0{Nc) due to trc, G^f is 0{g^N^) = 0{Nc), since gluon line in Fig.2 is 
equivalent to double fundamental line, yielding two color traces in Fig.2. 

From ( |45| ) one can see that the basic element defining the amplitude of the 
mixing of meson and hybrid is the dimensionless ratio 



2fig{n)\Mi''^ - Mm 



- Ar (50) 

The rest of this section is devoted to the calculation of the matrix element vj^l 
using realistic wave functions for the meson and hybrid. 

To make estimates of V/v/n, we first remark that to the lowest approximation 
in spin splittings both (fM and ^ipn are proportional to unit matrices in Lorentz 
indices, and moreover the gluon Green's function is proportional to 6^^^, 
where each bound gluon acquires its mass due to the attached string (similarly 
to the W^, acquiring mass due to attached Higgs condensate) and hence the 
sum over fii/ is the sum over 3 polarizations. 

As a result V^m^KiM ~ ^f^i^ ^^^^(0'^) does not depend on fi in the same 
lowest approximation (when gluon spin contribution to the hybrid mass is ne- 
glected; both for mesons and hybrids the spin splitting is less or about 10% of 
the total energy and our estimates will have this accuracy). 
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The meson wave function (pM{r) for linear growing potential is Airy function, 
for our purposes it is enough to use a Gaussian with the proper radius rg, e.g. 



3/4 _3 

I e ^ 



Using tq as a variational parameter, one gets 



(r^) =rl = 0.725fm2. 



(51) 



(52) 



For the hybrid wave function one can use the eigenfunction of the Hamilto- 
nian ( PI3| ) in the lowest hyperspherical approximation, where for an estimate 
we approximate W{p) by the oscillator well C{p — pof, C = near its min- 
imal value (the accuracy of this procedure is around one percent for the cases 
considered, see [|2^])). 
Then one can write 



^(ri3,r32) = 



1 



det 



ri r2 rs 
R ^ r? 



1/2 



IT 



1/4 



, 2~^P^P0) 



(53) 



Here m is an arbitrary mass, used for dimensional reasons and disappearing 
from final answers, po is found from the minimum of W{p), ^'{p = Po) = 0, 
and w is expressed through W"{po), mcu'^ = W" (po), so that 



UJ 



= (1.358)2/3^3^2/3 (i^l±J^ 



1/3 



2 //:(£ + 1; 



1/3 



P\Pz 



1/6 



m \ 1.358(7 



det 



pi + P3/ 

p = pi+ P2 + P3. 



ri r2 r3 
R ^ r/ 

The determinant ( |5B|) appears due to the change of variables, as follows 



(54) 

(55) 
(56) 



I (iri(ir2(ir3|^|2 = J dRd^dr]\ilj\^\ det {)\ = J dRp^dpdUQl'^fldetOl. (57) 

Finally / dQg = Qq = tt^ and we need expression of p through ri3,r23 etc., 
namely 



i<j '^P 



(58) 
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In Vmu enters ^„ at ri3 = 0, r23 = r2i = r, which yields 

2m/i 



/ = r^ ^''f + (59) 



Insertion of (|5T|)-(|5^) in yields 



VMn cg'^^ (60) 



where following notations are used: 

J = t'dte'''-^^'-'^\ (61) 

7r4(/:(/:+l))7/835/8l, ^3 ; [^) ' ^^^^ 

2a;r^/ii(^i + ^3) I Sflm 
a = — , b = po. 2 ■ -. 63 

In ( I6DI ) fig is the gluon effective energy, where gluon belongs to the hybrid, 
so that fig = fi^. 

We can do now estimates of Vmou both for light and heavy quarks. Using 
[PPI , |I| for massless quarks and for a = 0.18 GeV^ one has fi2 = fJ^i = 0.312 
GeV, /i3 = 0.442 GeV. u = 0.956 GeV, a = 2.54, b = 0.813, and hence 

= g ■ 0.2 J < g ■ 0.08 GeV. (64) 

For heavy quarks, e.g., for mi = m2 = 4.8 GeV, one has tq = 0.356 fm, 
/ii ~ mi = 4.8 GeV, fi^ = 0.767 GeV, and u = 0.5 GeV. Inserting all these 
values in (^D]) one obtains again for heavy meson-heavy hybrid mixing 

= g ■ 0.08 GeV. (65) 

This estimate is close to the calculations made in []5|, 0]. Hence one obtains a 
small mixing parameter A„ in ( |5tl|) when mass difference between meson and 
hybrid is large, A ~ 1 GeV, namely 

A„<0.l5 (66) 

while for close values of meson and hybrid masses, AM ~ 0.1 GeV, the mixing 
can be large and one should use the many-channel Hamiltonian, which can be 
obtained from (\1X\). 
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Figure 3: Glueball Green's function describing propagation of two gluons (wavy lines) from 
point X to point y, with insertion of quark-antiquark propagators (solid lines) from point u to 
point V. The cross-hatched region refers to the adjoint string world sheet, while hatched region 
between quark lines has the same meaning as in Fig. 1. 

5 Mixing of glueball with hybrid and mesons 

Here we consider a two-gluon glueball state, and using the same methods as 
discussed in Section 4 and Sections 1,2 we write the amplitude corresponding 
to the diagram of Fig. 3: 



G(f ^'^''^'(x,?/) = /tre{r(°-*)G^,(x,w)nr«(7,5'(7i,^)7,5'(z;,w))x 

xt^G',,(^,2/)r(i")G^v(x,2/)}. (67) 

Representing each gluon line at large N(. as a double fundamental line, one 
gets two closed color contours in Fig. 3 and as a result the factor A^^, hence 
the amphtude of Fig. 3 is 0{g'^N^) ~ 0{Nc) as well as the leading two-gluon 
diagram G^^J which means that mixing here also appears in the leading order 
of the 1/Nc expansion. 

Using expressions for the Green's functions (^, one obtains that in- 
sertion of the quark loop in Fig. 3 contributes to the glueball Green's function 
the term 

^ n 2iig{M\j' - Mg)[Mg- M\j') 

(68) 

Here (f^'^^fig is the effective (constituent) mass of the gluon in the initial and 
final glueball (which is calculated through dadj), Mq is the glueball mass, and 
other notations are the same as in the previous section. 
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Figure 4: The amplitude (Green's function) for transition of a glueball into a meson. Symbolic 
notations are the same as in Figs.l and 3. 



Similarly to the case of meson-hybrid mixing we introduce the mixing pa- 
rameter 

= A^", (69) 



l2,,,{G)\M'ff^ - Mg\ '• ' 

where Vmu is the same as in meson-hybrid case, (|55j) with the replacement of 
meson to the glueball Green's function 

g I (/?G(r)V^+(ri2 = 0,ri3 = = y)Sy. 



Mn 



(70) 



To estimate Vun one can use the same with following replacements: 

(i)/i^ is the effective mass (energy) of the gluon in the glueball, for lowest 
glueball states = 0.53 GeV (L = = 0), 0.69 GeV (L = 1, = 0) [M H; 

Insertion of the values of ^u-i, ;U-2, /U3, cj, into a, 6 yields a = 1.17a, h = 
0.9236, hence all estimates for Vun made in (^, ( |55| ) hold also for Vun, and 



the glueball-hybrid mixing is the same as meson-hybrid one, provided the radius 
ro of glueball is the same as that of meson. 

We come now to the amplitude of Fig. 4, which gives the amplitude for 
glueball-meson mixing. Analysis similar to the previous one yields for this 
amplitude the following ratio: 



\GM _ 
A„ — 



nM 



1 



in) 



_ \GH\HM 

— K K ■ 



71: 



Insertion of our results leads to the following order of magnitude 

estimate 

^2(0.08 GeV)2 



A 



GM 



{AMY 



(72) 
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For AM ~ 1 GeV one obtains A^^^ ~ 1/16, and the probability of admixture 
of the glueball to a meson Pgm = l-^n^'^^P is less than a percent. This should be 
compared to the probability of the hybrid in the same meson 

Phm = lA^T = iVc/ (^^1^)' - ' ~ 5 - 10%. (73) 

All these estimates hold for unsuppressed glueball-hybrid-meson transitions, 
when spin-flip amplitudes do not enter, otherwise one may expect additional 
one order-of-magnitude suppression of probabilities. 



6 Discussion and conclusions 



In Sections 4 and 5 we have obtained reduced matrix elements for meson-hybrid, 
hybrid-glueball, and meson-glueball transitions. The full matrix elements in- 
volve expressions of the following type (see (^) and (|68|)): 



A-r = /{$(f')J,7.*£'dr3, (74) 

(H) 

where contains all components of the hybrid wave function. 

Since the Hamiltonian ([TD|) and the eibein technic (|TT|) separate only posi- 
tive values of fiq, /i^, and /ig, the wave functions and refer to those 
positive values and we disregard the negative energy components (negative /i- 
components). 

Moreover, spin interactions in mesons, hybrids, and glueballs are treated 
perturbatively in our formalism, and comparison to lattice data for the same 
input shows that both approximations of positive fii and perturbative spin 
splittings are rather good for the cases considered - see for Tables with 
comparison. These approximations are not good for pion and kaon, where both 
spin and negative energy states play an important role, but we shall not consider 
these mesons, or hybrids with similar properties. 

Therefore our parameters (^D|) refer only to such states which are 

well described by the Hamiltonian (|TD|) with condition (\n\). In this case the 
matrix element of 7^ is simplified since the fi dependence in the hybrid decouples 
from the spin (Lorentz) indices and the simplifications made after ([5D|) are valid, 
when one neglects possible Clebsch-Gordon coefficients, which are of the order 
of unity. 

In this way one obtains an estimate for the "superallowed" transitions meson- 
hybrid-glueball, not involving spi-flip or negative-sate components. 
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A comparison with the results of [j^, shows a good agreement within a 
factor of two, which tells that the missing in our case factors are of the order 
of unity. However, in 1^,0] as well as in our case only wave functions were 
considered without negative energy components (negative p-spins). 

To have more comparisons, one may look at the lattice calculation of MME, 
which have been done extensively for the glueball 0++- scalar meson case 



The authors []2^ arrive at the following result of careful lattice studies: 



|/o(1710) >= 0.859|£/ > +0.302|ss > +0.413|nn >, 

|/o(1500) >= -0.128|i/ > +0.908|ss > -0.399|nn >, 

|/o(1390) >= -0A9b\g > +0.290|ss > +0.819|nn > . (75) 

^From ([75| ) it is clear that a strong mixing occurs between states in the 
region 1.4-1.7 GeV. However our results for MME between glueball and meson 
states always implied an appearance of a hybrid state as an intermediate state 
between glueball and meson. From (|7I|) it is clear that the strong glueball- 
meson mixing is possible only if the intermediate hybrid has the mass in the 
same interval 1.4-1.7 GeV introducing in ( |72| ) AM ~ 0.2 — 0.3 GeV, one obtains 
^GM ~ 1^ as in the lattice calculations []2§] for the 0++ states. 

Hence one could look for the 0"''+ hybrid state in the discussed mass range. 
Indeed, analytic calculations of lowest hybrid states in confirm the possi- 



bility of the 0++ state around 1.7-1.8 GeV. this problem clearly calls for further 
theoretical and experimental investigation. 

The author is grateful to K.A.Ter-Martirosyan for many stimulating discus- 
sions, suggestions and criticism of the first drafts of the paper. This work was 
partially supported by the RFFI grants 00-02-17836 and 00-15-96786. 

Appendix 
Meson wave-functions with definite 

In (PT|) the wave function of the meson was given in the so-called 4x4 
representation, where the total angular momentum J and parity P were not 
specified. 

In this Appendix some relations and notations are given for meson wave 
functions with given J^, based on last reference [j2^ . 
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The states of the quark may be classified according to the so-called p-spin 
and usual spin states, where /o-spin states, p = ±1, can be taken as eigenvalues 
of each of 70*^ (i = 1,2 refer to quark 1 and quark 2, which may be an antiquark), 
or as eigenvalues of A^, where 

A'^(p) = p(7P + ^)+7oVp^ + ^^ _ 
In both cases one can create four p-spin eigenstates (e.g., as eigenstates of 

(1) (2)n 

To To ): 

| + + >-|--> ! + +>+!--> 

a 



\/2 ' V2 

Using one can form 8 states of unnatural parity for qq, P = {—)'^~^^ 



B = {'J}, 'rj, V}, V^, ^(J-l)}, ^(J-in, '(J+l)}, 'U+l^}, (A.3) 
and 8 states of natural parity P = {—Y: 



B^ = {'j}. 'Ji V}, v-m, 'u+m, %j+m- (A.4) 

Dynamical equations for mesons connect all components belonging to the 
same . 

The 16, component wave function (PT|), f^^^'^, can be decomposed into 4x4 
qq basis, introduced in the last reference [)2^, as it is shown in the Table below. 

A similar expansion can be derived for the hybrid states, which include in 
general 16 x 3 components. In the estimates of the present paper it is assumed 
that the positive energy (positive p-spin) component gives the dominant contri- 
bution to the mixing matrix element, and that spin interaction can be treated 
perturbatively. Then the estimates (^, (^ are modified due to kinematical 
coefficients (like those in the Table) by factors of the order of unity. 
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Table 1: 16-component vector gg-basis states and corresponding 4x4 matrix gg-basis states. 
(The angular dependence of the wave functions of the gg-triplet states must determine to 
which gg-triplet state they correspond. For a correct normalization, all matrix states should be 
multiplied by 1/2; furthermore, matrix states proportional to p or p need an extra factor 1/M 
and l|p|, respectively. The correspondence is only valid in the c.m.s. and po = is assumed.) 
Here P and p refer to the total and relative momentum of gg in the meson respectively. 
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